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Abstract 

We derive an exact Green's function of the diffusion equation for a pair of 
spherical interacting particles in 2D subject to a back-reaction boundary condi- 
tion. 

1 Introduction 

Green's functions (GF) for an isolated pair of molecules play an important role in 
the theory of diffusion-influenced reactions for several reasons. First, the solution 
for any other initial distribution can be calculated when this GF is known [21 
[3l [7]. Second, the GF can be used to derive important other quantitities, for 
instance the survival probability and the time-dependent rate coefhcient [2 [3j [7] . 
Third, the GF may be used in particle-based stochastic simulation algorithms to 
enhance the efficiency of naive Brownian dynamics simulations, cp for instance 
[5]. Finally, the knowledge of an exact analytic expression permits to validate 
newly devised stochastic simulation algorithms [S]. 

Exact analytic expressions have been derived for the GF of an isolated pair 
that can undergo a reversible reaction in ID and 3D [3J [31 [7] . To the best of 
our knowledge, no exact expression has been presented so far for the 2D case. 
However, in particular with regard to cell biological applications, an accurate 
theoretical treatment of two-dimensional diffusion-influenced reactions provides 
the basis for a better understanding of important processes like receptor cluster- 
ing on cell membranes [1]. 

To derive the GF in 2D we consider an isolated pair of two spherical particles 
A and B with diffusion constants Da and Db, respectively. The particles may 
associate when their separation equals the "encounter distance" a to form a 
bound molecule AB. When bound, the molecules may dissociate again to form 
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an unbound pair A + B. Such a system may be described as the diffusion of a 
point-like particle with diffusion constant D = Da + Db around a static sphere 
with radius a. In this picture, reactions are introduced by imposing boundary 
conditions at the sphere's surface. The association reaction is described by the 
radiation boundary condition that is characterized by an intrinsic association 
constant Kq. The radiation boundary condition is used to describe irreversible 
association. To take into account reversible reactions, i.e. dissociations also, 
the radiation boundary condition has to be generalized to the back-reaction 
boundary condition that involves an additional intrinsic dissociation constant 

We consider the probability density function g{r,t\rQ) for the probability to 
find the particles at a distance equal to r at time t given that the distance was 
inititially at time i = 0. The time evolution of g(r, t|ro) is governed by the 
2D diffusion equation 

^^g{r,t\r^)^D{^ + ~)g{r,t\r^), r>a (1.1) 

The diffusion equation has to be completed by specifying boundary conditions. 
Together with the following initial 

2TTrg{r,to\ro,to)=6{r-ro) (1.2) 

and boundary condition 

(7(r-^«),t|ro,to) -0 (1.3) 

equation (jl.ip (if a = 0) is equivalent to the free-space diffusion equation in 2D 
with the familiar solution 

<7free(r,i|ro) = ^ ('^'+'-0 )/4^*/o ( ^ ) (1.4) 

also known as the free-space Green's function. 

However, as described previously, the PDF we are interested in is only defined 
for r > a > and one has to impose a boundary condition for r = a specifying 
the behavior at the encounter distance. We will discuss the following case [5] [7] : 

d 

2TTaD—g{r,t\ro)\r=a = i^agir = a,t\ro) - Kd[l - S{t\ro)]. (1.5) 
or 

Here S(t\rQ) denotes the survival probability that a pair of molecules with initial 
distance ro survives by time t 

Sit\ro)^l- f 27TaD-^g{r,t'\ro)\r=adt' (1.6) 

Following [5], we make the following ansatz for the Laplace transform of the 
Green's function that satisfy the backreaction (BR) boundary condition 

9{r, q\ro) = 5froo(r, giro) + gsnir, q\ro). (1.7) 

Here 

ff.ee(.,g|ro) = ^| /o(<zr)i^o(<Z.o) r<ro ^''^^ 
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is the Laplace transform of the free-space Green's function (jl.4|) . The variable q 
is defined by 

9 := V^, (1-9) 

where p denotes the Laplace domain variable. The part ^br that takes into 
account the boundary condition is a solution of the Laplace transformed 2D 
diffusion equation 

, 2 H 1 q 9BR = 0. (1.10) 

ar^ r dr 

The general solution to p.lOp is AKo{qr) + BIo{qr), where Iq{x), Kq{x) refer 
to the modified Bessel functions of first and second kind, respectively, and of 
order zero Because we require lima;_>.oo 3br 0, and limx-^oo Ia{x) oo, 
the coefficient B has to vanish and hence, 

5BR(r, q\ro) = A{q, ro)Ko{qr). (1.11) 

A{q,ro) is determined by the requirement that the complete Green's function 
(|1.7I) satisfies the Laplace transformed backreaction boundary condition, cp. 

d d 

^5br(?', 9 ko) I r=a = hgBR{r,q\ro)\r=a - HdP^^—gBR{r,q\ro)\r=a, (1-12) 
where we have defined h := Using (fTTl) . (|LS)) . dTTTTl) . (|TTT^ and 

/^(x) = /i(x), (1.13) 
K^{x) = -K,{x) (1.14) 

and defining kd :— jj, we obtain 

~ / I N 1 {q^ + KD)Ii{qa) - hqIo{qa) 

9BRir,q\ro) = ^ + .^)i,,(g,) + /,,;^„(,,) Koi,ro)Koiqr). (1.15) 

The inversion theorem for the Laplace transformation can be applied to find 
the corresponding expression in the time domain 

9BRir:t\ro) = — eP* gBRir, q\ro)dp. (1.16) 

To calculate the Bromwich contour integral (|1.16l) we first note that ^br has a 
branch point at p = 0. Therefore, we use the contour of Figure 1 with a branch 
cut along the negative real axis, cp. [5]. Furthermore, because the integrand has 
no poles within and on the contour [6l and because the contribution from the 
small circle around the origin vanishes, we obtain 

e^* .9BR(r-, q\ro)dp = / e^* ffBR(r, q\ro)dp + 

J 7 — 100 

+ / eP'gBR{r,q\ro)dp+ / gBR{r,q\ro)dp. (1.17) 

JC2 JCi 



Thus, it remains to calculate the integrals , 
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We now choose 

p^Dx^e""^ (1.18) 

and use [5] 

/„(xe±^'/2) = e±""/V„(x), (1.19) 
i^„(xe±"/2) ^ ±i7rzeT""/2hJ„(x)±iy„(x)]. (1.20) 

Jn{x),Yn{x) denote the Bessel functions of first and second kind, respectively 
[T]. It fohows that 

n ■ nOO 

eP'gBR{r,q\rf,)dp^- e-""-"' H^^\xr)H^^\xro)F{x)xdx. (1.21) 
(2) 

Here Hn {x) := J„(a;) — iYn{x) denotes the Bessel function of third kind (also 
referred to as Hankel function) 1 and we have defined 



p/ N _ a{x)[aix) +iPix)] 



and 



a{x) :— {—x^ + K£j)Ji{xa) — hxJo{xa) (1-23) 
P{x) := {-x"^ + KD)Yi{xa) - hxYoixa) (1.24) 

To evaluate the integral along the contour C4 we choose p = Dx^e^™ and 
after an analogous calculation one finds that 

e^* gBRir, q\ro)dp = - eP* gsRir, q\ro)dp (1.25) 

C2 JCi 



where J^^ means complex conjugation. Thus, one arrives at 

9BK{r,t\ro) = TT- j 9BRir,q\ro)dp (1.26) 

= -^^{l eP'~gBR{r,q\ro)dp^ (1.27) 

_ I r _a,2 Mx)[a{xMx) + l3{x)U{x)] 

2 Jo a{xY+P{xY ^'-'^^ 



where we have defined 



n{x) := Jo{xr)Jo{xro) - Yoixr)Yo{xro) (1.29) 
n(a;) := Yo{xr)Jo{xro) + Jo{xr)Yo{xro). (1.30) 



Next, we use the fact that the free-space Green's function may be written as 



1 /"^ 9 

gfrco(f,t|ro) = — / e^"'''*jQ{xr)Jo{xro)xdx (1.31) 
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to arrive at the exact Green's function in the time domain 



/ 

Jo 



■ oo 



g{r,t\ro) = — 

2,71 



Dx^t 



To{xr)To{xro)xdx 



(1.32) 



e 



with 



To{xr) 



Jo{rx)[{x'^ — K, D)y^i {xa)-\-hxYo{xa)]—Yo {rx) [{x"^ — K]j)Ji {xa)-\-hxjQ{xa)] 
{[{x'^-KD)Ji{xa)-\-hxJo{xa)]'^-\-[{x'^-KD)Yi{xa)-\-hxYo{xa)]^y/'^ 



(1.33) 



Note that the hmit — > one recovers the known Green's function ^ for the 
irreversible case with radiation boundary condition. 
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